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Abstract 

We show that the intracavity Kerr nonlinear coupler is a potential source of bright continuous 
variable entangled light beams which are tunable and spatially separated. This system may be 
realised with integrated optics and thus provides a potentially rugged and stable source of bright 
entangled beams. 

PACS numbers: 42.50.-p,42.50.Dv,42.65.-k 
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I. INTRODUCTION 



The term nonlinear coupler was given to a system of two coupled waveguides without 
an optical cavity by Perina et al. [lj]. Generically, the device consists of two parallel optical 
waveguides which are coupled by an evanescent overlap of the guided modes. These may 
operate with both modes propagating in the same direction, or in opposite directions. The 
quantum statistical properties of this device when the nonlinearity is of the x^ type have 
been theoretically investigated , among the predictions being energy transfer between 



the waveguides |4( and the generation of correlated squeezing 



The production of entan- 



gled states of the electromagnetic field with an intracavity Kerr nonlinear coupler in which 
one of the modes is externally pumped was analysed by Leohski and Miranowicz |6| . In the 
travelling wave configuration, the generation of entanglement via coupled parametric down- 
conversion has been analysed by Herec et al. [7\. Olsen and Drummond [8] and Olivier and 
Olsen have analysed coupled intracavity downconversion in terms of continuous variable 
entanglement in both the above and below threshold regimes. The system we consider here 
consists of two coupled nonlinear materials with x nonlinearity, operating inside a pumped 
Fabry- Perot cavity. The details of the coupling are as in the work of Bache et al. [10| , who 
analysed a scheme using coupled second harmonic generation in terms of intensity correla- 
tions between the modes. The coupling is realised by evanescent overlaps of the intracavity 
modes inside the nonlinear medium, which can be either a single Kerr crystal pumped by 
two spatially separated lasers, or two Kerr waveguides. This t ype of coupling has already 
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been investigated both theoretically and experimentally 

In this paper we will calculate the phase- dependent correlations between the two outputs 
of the cavity which demonstrate that the system is a source of continuous variable entangled 
states and states which exhibit the Einstein-Podolsky-Rosen (EPR) [1J| paradox. The spa- 
tial separation of the entangled output modes means that they do not have to be separated 
by optical devices before measurements can be made, along with the unavoidable losses 
which would result from this procedure. The correlations are tunable by controlling some of 
the operational degrees of freedom of the system, including the evanescent couplings between 
the two modes, the input power and the cavity detunings. The entanglement is also between 
outputs which can be of macroscopic intensities, which may be an advantage compared to 
the nondegenerate optical parametric oscillator (NDOPO) operating above threshold. Vil- 
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lar et al. |15( recently demonstrated bright entangled outputs with such a device, but were 
not able to make an unambiguous demonstration of the EPR paradox due to losses. The 
scheme we analyse here may provide a more stable route to a demonstration with spatially 
separated bright beams. 



II. HAMILTON! AN AND EQUATIONS OF MOTION 

We consider two evanescently coupled x materials inside a pumped Fabry-Perot, which 
could be experimentally realised with, for example, a dual core fiber with dielectric mirrors 
at each end The interaction Hamiltonian in the rotating-wave approximation can be 
written as 

7~~£int 7~~tpurnp 7~~tdamp 'Hherr ^-couple} (1) 

where the pumping Hamiltonian is 

2 

ftpump = ihy^ j fej-qt - €*dj) , (2) 

3=1 

the damping Hamiltonian is 

2 

Hdamp = ( r i«i + > ( 3 ) 

3=1 

the Kerr Hamiltonian is 

2 

H kerr = H ^ Xj^j 2 , (4) 
3=1 

and the coupling Hamiltonian is 

Hcoupie = ft-J {aid\ + a\a<^ . (5) 

In the above, the dj are the annihilation operators for the intracavity modes, the 6j represent 
the classical pumping terms, the Tj are bath operators, the Xj represent the Kerr nonlin- 
earities, and J represents the evanescent coupling. The form of the latter is as described in 
Bache et al. |lf| . 

As we do not expect to be able to solve the Heisenberg equations of motion, we will derive 
stochastic differential equations in the positive-P representation [17], after making the usual 
zero-temperature Markov and Born approximations for the cavity damping Proceeding 
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[19], 



via the usual methods [19], and making the correspondences between the operators dj, 
and the c-numbers (Xj, ad", we find the following Ito equations, 

dai 
~dt 
daf 



ei - (71 + iAi)ai - 2ixiafal + iJa 2 + \ —2ixi<x\ Vi 



dt 

da 2 



dt 
da 



e i _ (7i _ i^i)af + 2iXi a f 2 «i - iJol% + y 2ixiaf 2 r} 2 , 
£2 - (72 + iA 2 )a 2 - 2ix 2 a\ «2 + iJ a i + v/-2ix 2 al r] 3 , 



- ~~ e 2 - (72 - iA 2 )aJ" + 2ix2at 2 «2 -iJaf + J 2%x 2 a\ 2 »?4, (6) 



where the r)j are real Gaussian noise terms with the correlations T]j(t) = and T)j(t)r)k(t') = 
5jkS(t — t'). We now have a set of equations which are a coupled version of those used 
by Drummond and Walls to describe a model of optical bistability with an intracavity Kerr 
medium j^Q]. In principle these equations may be solved to calculate any quantity which can 
be written as a time-normally ordered operator product. However, rather than immediately 
solving the full equations, we will linearise them about their classical steady-state solutions 
and examine the stability of the system. This is equivalent to treating the system as an 
Ornstein-Uhlenbeck process j^j] and allows for simple determination of the output spectra 
in parameter regions where the process is valid. 



III. LINEARISATION AND STABILITY 



The classical equations are found by dropping the noise terms from Eq. |U1 and are 

ei - (71 + iAi)ai - 2ixi\cei\ a\ + iJa 2 , 



dt 

—^T = e 2 - (j 2 + iA 2 )a 2 - 2ix 2 \a 2 \ 2 a 2 + iJa t . (7) 

To begin with, we will consider that all the parameters are symmetric, that is e\ = e 2 = e, 
7i = 72 = 7, Ai = A 2 = A and Xi = X2 = X- We now write equations of motion for 
a p (= «i + a 2 ) and a m (= a\ — a 2 ) as 

da 



dt 
~df 



v - = 2e- [j + i(A- J + 2 X \a\ 2 )] a p , 

-[ 1 + t(A + J + 2 X \a\ 2 )]a m , (8) 
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where we have used the symmetry of the system to set = \ot2\ = \ a \ ■ The steady-state 
solutions of the above equations give us af% = 0, so that af = af = a = a s p s /2 and 

a = [ 7 + i(A- J + 2xH 2 )]" (9) 
Multiplying this equation by its complex conjugate and setting e G 9ft, we find a cubic 



equation for I = \a 



2 



4 X 2 / 3 + 4(A - J) X I 2 + [ 7 2 + (A - J) 2 ] / - e 2 = 0. (10) 



A condition for bistability is that the quadratic equation found by differentiating Eq. El 
with respect to / has two positive real roots. The quadratic is 

12 X 2 / 2 + 8(A - J) X I + 7 2 + (A - J) 2 = 0, (11) 

with roots 



r± = — 
6x 



2(J-A)± v/(J-A) 2 - 3 7 2 . (12) 
The condition then is that (J — A) 2 > 37 2 , so that bistability is possible provided that 

J>A + v / 37- (13) 

We note here that this condition is necessary but not sufficient, as it does not tell us at 
which values of the pumping the bistability may be found, for example. We also note that 
for the case where A = J, similar to that which was used to simplify the analysis for coupled 
downconverters j^, 0] , there is no bistability. It is obvious that for the cavity on resonance 
(i.e A = 0), bistability is possible for any value of the coupling greater than a/37, as shown 
in Fig. and that for A < — V07, bistability is possible for any value of the coupling. 
We will not investigate this bistability further at this time except to remark that a fully 
quantum analysis of an uncoupled oscillator showed that quantum fluctuations removed the 
bistability which was seen in a semi-classical analysis (^(J. 

We will now return to the steady-state solutions themselves. Eq. El m &y be solved ana- 
lytically to give solutions which are rather complicated, so that we will proceed numerically 
in the general case. However, for A = J, the real solution simplifies to 

2/3 

(14) 



-3 2 / 3 7 2 + 3 1 / 3 (9 X e 2 + ^3(76 + 27e 4 x 2 

1 = 1 \ 1/3 

6 X (9%e 2 + V3( 7 6 + 27e 4 X 2 )) 
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FIG. 1: The intensity as a function of pump intensity, for % = 10 , 7 = 1 and A = 0, and two 
different values of the coupling. All quantities shown in this and subsequent figures are dimension- 
less. 

which may then be substituted into Eq. |U] to give an expression for the mode amplitudes. 
The values thus found agree with the steady-state mode amplitudes found by numerical 
integration of the equations of motion. 

Now that we have found steady-state solutions, we may linearise the equations of motion 
around these by setting ctj = a + 5a j and expanding to first order in the 8a j. This gives 
the equations 

— Sax = —(7 + iA)5ai — 2i\ (a 2 5 at + 2|a| 2 <5ai) + iJ5a 2 + \/—2ixa 2 m 



_ 5 a + = -(7 - iA)6a% + 2%x (a* 2 5a 2 + 2\a\ 2 5a 2 h ) - iJ5a\ + ^2%xa* 2 774, (15) 



(7 - %A)8a\ + 2ix (a* 2 5a 1 + 2\a\ 2 5af) - i.J5a\ + ^/2%xa* 2 r} 2 , 



(7 + iA)5a 2 — 2%x (a 2 5a 2 + 2\a\ 2 5a 2 ) + iJ5a\ + a/— 2%xa 2 773 




which may be written in matrix form as 



-r5a = — Aba + B, 



(16) 



where 5a = [Sai, Saf, 5a 2 , Sa^] 7 , the drift matrix is 
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7 + i(A + 4x|«| 2 ) 2ixa 2 -iJ 

-2ixa* 2 7-i(A + 4xM 2 ) 

-iJ 7 + i(A + 4x|a| 2 ) 







U 



-2ixa 



* 2 





U 
2i X a 2 
7-2(A + 4x|«| 2 ) 



(17) 

and BB T is a diagonal matrix with the vector [— 2i\a 2 , 2ixct* 2 , — 2ix<y 2 , 2i\a* 2 } on the 
diagonal. The validity of this process depends on the eigenvalues of the matrix A having no 
negative real part, which we will verify numerically in all reported results. These eigenvalues 
may be found analytically as 



Ai, 2 = 7 ± v^xM 2 [3xl«| 2 + 2(A - J)} - ( J + A)2, 
A 3i4 = 7 ± v / 4x|«| 2 [2(J-A)-3 X |«| 2 ]-(J-A) 2 . 



(18) 



For A = J, which we will use in the results presented below, the linearisation process is not 
valid for |a| 4 > (7 2 + 4J 2 )/12% 2 . In this region, we would need to use stochastic integration 
of the full equations (jOJ to obtain results. 

Writing the equation in the above form (fTT)|) leads to a particularly simple expression for 
the intracavity spectral correlation matrix, 



S(u) = (A + iuiy 1 BB T (A T - iwl) 



-i 



(19) 



after which we use the standard input-output relations [22j to relate these to quantities 
which may be measured outside the cavity. The matrix S(u) contains all the information 
necessary to calculate the correlation functions which we will give in section IIVI 



IV. ENTANGLEMENT CRITERIA 



Entanglement is a property of quantum mechanics which is related to the inseparability 
of the combined density matrix of a system into density matrices for its subsystems. In the 
present situation, we are interested in continuous variable bipartite entanglement between 
the modes 1 and 2. The presence of this entanglement can be shown in a number of ways. 
Before we proceed with the definition of the criteria we will use here, we will define the 
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optical quadratures which will be used, as their actual definition affects the form of the 
inequalities. We define the quadrature operators at the phase angle 6 as 

X* = a,e- w + a]e w , (20) 



so that [Xj, xl +n ^ 2 ] = 2i5jk and therefore the Heisenberg uncertainty principle (HUP) 
requires V(Xj)V(Xl +7T ^ 2 ) > 15jk- In the interests of notational simplicity, we will label the 
quadrature X 6 k +7T ^ 2 as Y k e . We also note here that our results will be presented in terms of 
output spectral correlations, as these are the quantities that are experimentally measurable, 
and that the same inequalities as presented below will hold for these. We will also present 
the correlation functions in terms of intracavity variances, again for reasons of notational 
simplicity. 



2J, who 



The first of the entanglement measures is due to Duan et al. 23] and also Simon 
developed inseparability criteria which are necessary and sufficient for Gaussian states, and 
sufficient in general. These criteria have recently been shown to be special cases of an infinite 
series of inequalities based on the non-negativit y o f determinants of matrices constructed 



from certain combinations of operator moments 25[. In the general case, we may define the 
quadrature operators similarly to Duan as 

*l = \b\Xl ± ^Xl 

n = \b\Y*±±fi, (2i) 

where b is an arbitrary non-zero real number. It may be shown that, for separable states, 

V(Xi) + V(Y°) > 2 (V + i) , (22) 

with any violation of this inequality therefore demonstrating the presence of bipartite en- 
tanglement. In what follows, we will choose b = 1 so that the lower bound of the inequality 
is 4. While this is not the optimal choice for the general case, it is sufficient for purposes 
of comparison with the other measures we will use to categorise our system. In Fig. |21 we 
give the predictions for the values of the spectral output correlation S out (X e _) + S out (Y+) for 
different values of the interaction. We see that, as the interaction strength increases, the 
maximum violation moves away from zero frequency and is found at different quadrature 
angles, with the spectra becoming bifurcated. This is different from the cases with coupled 
interactions, where setting the detunings equal in strength to the evanescent couplings 



meant that the optimal correlations were found for zero quadrature angle 0,0]. It happens 
because the x^ interaction itself also changes the quadrature angles at which the optimal 
correlations may be measured J26], acting in some sense as an intensity dependent detun- 
ing [2^. In principle it may be possible to dynamically control the detuning, coupling and 
intensities so that the largest violations of the inequalities were found for 8 = 0, although we 
will not investigate this further because it is relatively simple to control the phase of local 
oscillators to choose the appropriate angles. We have shown the results for the quadrature 
angles at which the violation of the inequalities is maximised for different values of the non- 
linear interactions. The existence of the maximal violations at finite frequencies could be 
a real operational advantage, due to the excess noise which is found experimentally at low 
frequencies. 
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FIG. 2: Violation of the Duan inequalities for e = 10 3 ,7 = 1, J = A = 10. The solid line is for 
X = 10~ 5 , at a quadrature angle of 6 = 80°, the dashed line is for x = 10 -6 , at an angle of 9 = 122° 
and the dash-dotted line is for x = 10 -7 , at an angle of 9 = 14°. The inequality is violated and 
entanglement is demonstrated for any value less than 4. 



An alternative to indicate the existence of bipartite entanglement is a demonstration of 
the EPR paradox, which, as proven by Reid is sufficient to demonstrate inseparability 



of the system density matrix. The EPR paradox itself was introduced in an attempt to show 
that quantum mechanics was not a complete and locally realistic theory Schrodinger 
replied by introducing the concept of entangled states which were not compatible with 
classical ideas such as local realism j^. In 1989 Reid proposed a physical test of the EPR 
paradox using optical quadrature amplitudes jsjj, which are mathematically identical to 
the position and momentum originally considered by EPR. Reid later expanded on this 
work, demonstrating that the satisfaction of her two-mode 1989 criterion, and hence a 
demonstration of EPR correlations, always implies bipartite quantum entanglement . It 
was also shown by Tan jsjj that the existence of two orthogonal quadratures, the product 
of whose variances violates the limits set by the Heisenberg uncertainty principle, provides 
evidence of entanglement. Tan demonstrated this in the context of teleportation, with the 
outputs from a nondegenerate OPA mixed on a beamsplitter. 

To examine the utility of the system for the production of states which exhibit the EPR 
paradox, we will follow the approach of Reid (^J, as outlined in Dechoum et al. We 
assume that a measurement of the Xf quadrature, for example, will allow us to infer, with 
some error, the value of the X\ quadrature, and similarly for the Y® quadratures. This 
allows us to make linear estimates of the quadrature variances, which are then minimised 
to give the inferred variances, 

v(xlx°)X 

V*'(X') = V(X°) - 1 J , (23) 

V(Yf,Y 2 d )] 2 

V*»'{Yf) = V<Yf) - 1 J , (24) 

where V(A,B) = (AB) — (A)(B). The inferred variances for the j = 2 quadratures are 
simply found by swapping the indices 1 and 2. As the X e - and Y® operators do not com- 
mute, the products of the actual variances obey a Heisenberg uncertainty relation, with 
V (Xj)V iXj) — 1- Hence we find a demonstration of the EPR paradox whenever 

V inf (X<j)V inf (Yf) < 1. (25) 

In Fig. El we show results for the product of the inferred spectral output variances, indicating 
a clear demonstration of the EPR paradox and hence entanglement between the modes. The 
violations are found to be maximum for the same quadrature angles as found for the Duan 
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inequalities, with a qualitative similarity between the results. The small violations at higher 
frequencies, shown in Fig. El for the Duan inequalities, were not so apparent for the products 
of the inferred variances. This can be understood because the Duan criterion provides 
necessary and sufficient conditions for entanglement with a Gaussian system, while the EPR 
correlation provides a sufficiency criterion. Operationally, it would not be a problem as 
the major violations at lower frequencies are potentially more interesting, due to the higher 
intensities of the output fields closer to the cavity resonance frequency. 




(o (units of y) 

FIG. 3: Violation of the EPR inequality for e = 10 3 ,7 = 1, J = A = 10. The solid line is for 
X = 10 -5 , at a quadrature angle of 9 = 80°, the dashed line is for x = 10 -6 , at an angle of 9 = 122° 
and the dash-dotted line is for x = 10~ 7 , at an angle of 9 = 14°. 

The third measure which we will apply is the logarithmic negativity, proposed by Vidal 
and Werner as a com put able measure of entanglement, as opposed to others which can be 
difficult to calculate |33j. This measure has been used by Mancini to quantify the effects 
of feedback on continuous variable entanglement in a two-mode system [34]. We note here 
that this measure is defined for Gaussian states, to which we are also limited here due to 
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the linearisation process we have used. We first define the system covariance matrix as 



C 



C21 C-2 



(26) 



where 



and 



Defining 



IJ 



V{X°) V{X°,Yf) 
Vtff,X') V(Yf) 

'V(X°,X°) V(Xf,Yf 
V(Yf,X}) V(Y°,Y* 



£ = y (detCi - detCia) - y/(det C 2 - det C u ) 2 - det C 
the logarithmic negativity is then defined as 



(27) 



(28) 



(29) 



-io g2 e if e < 1 

otherwise. 



Any non-zero value of is then an indication that the two modes are entangled. An 
interesting feature of this measure is that it has no dependence on quadrature angle, at least 
with the system we are investigating here. The values shown in Fig.|3]are valid for all angles, 
even though the other two measures will not indicate entanglement for arbitrary angle. This 
shows that the logarithmic negativity is useful for demonstrating that continuous variable 
entanglement exists in a given system, but does not tell us at which quadrature angles the 
system may exhibit the necessary properties for uses such as teleportation. We note that 
the frequencies at which the logarithimic negativity is non-zero are the same as those at 
which the Duan inequalities are violated. 



V. CONCLUSIONS 



We have analysed the intracavity Kerr nonlinear coupler in terms of its entanglement 
properties and shown that it can produce macroscopically intense outputs which possess 
continuous variable entanglement. This entanglement has been demonstrated by the cal- 
culation of three different measures, the first two of which give maximum violations of the 
appropriate inequalities at particular quadrature angles, and the third of which indicates the 

12 




presence of entanglement in a manner which is not phase dependent. As this device could 
be constructed using integrated optics and only needs to be optimised for one frequency 
of operation, it could provide a rugged and stable source of bright entangled beams which 
could be of some utility for such technological uses as, for example, quantum communication 
channels and continuous variable teleportation. 
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